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Abstract 

The combined effect of viscosity and vorticity on the growth rate of the bubble associated 
with single mode Rayleigh -Taylor instability is investigated. It is shown that the effect of 
viscosity on the motion of the lighter fluid associated with vorticity accumulated inside the 
bubble due to mass ablation may be such as to reduce the net viscous drag on the bubble 
exerted by the upper heavier fluid as the former rises through it. 
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I INTRODUCTION 

Recently obtained simulation results^ regarding single mode Rayleigh Taylor instability(RTI) prob- 
lem indicates that the bubble developed at the two fluid interface is accelerated to a velocity quite 
above the classical model valued . This is ascribed to vorticity formation on the bubble-spike inter- 
face that is suggested^ to diminish the friction drag due to viscosity. The theoretical calculations 
and the consequent numerical results regarding Rayleigh Taylor instability in a viscous fluid- also 
indicates growth of vorticity inside the bubble in the neighbourhood of its tip, where the curvature 
is maximum. Recent experimental and simulation results^ show that the influence of viscosity on 
RTI may strongly suppress the growth rate. It is suggested that high pressure and strain rate condi- 
tions give rise to a phonon drag mechanism resulting in the lattice viscosity provided a solid state is 
maintained^. Moreover, it was also suggested by Betti and Sanz^ that at the spike, mass ablation 
induces a transverse component of velocity and thus gives rise to vorticity generation. It has been 
shown by the same authors that the instability growth rate is augmented by vorticity accumulation 
inside the bubble resulting from mass ablation. 

As the bubble rises through the denser fluid it encounters viscous resistance. In this paper we 
have incorporated the joint effect of viscosity and vorticity. It is shown how the vorticity induced 
contribution to bubble velocity as derived by Betti and Sanz^ is affected by the viscous drag exerted 
by the lighter fluid inside the bubble. It is found that the combined effect of vorticity and viscosity of 
the bubble fluid tends to diminish the friction drag exerted by the upper heavier fluid as the bubble 
penetrates through it. In this respect it is similar to the friction drag reduction effect due to vorticity 
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as suggested by the phenomenological model of Ramaprabhu et al.W . 

Section II contains the mathematical formulation as also the two fluid boundary conditions. The 
results and discussions are presented is section III. 



II MATHEMATICAL FORMULATION AND INTERFACIAL 
BOUNDARY CONDITION 

We are considering 2D single mode Rayleigh- Taylor instability problem. The y axis is taken vertically 
upward opposite to the direction of gravity and x axis taken along the horizontal plane of unperturbed 
two fluid interface plane. The equation to the perturbed surface of bubble which is assumed to be 
parabolic is 

y = r]( y x,t)=r] (t)+r ]2 ( y t)x 2 ; y>0 (1) 

with T]o(t) > and r] 2 (t) < (terms 0(x 2+t )(i > 1) are neglected as we are interested only in the 
motion very close to the tip of the bubble \x\ <C 1)®. 

The kinematic boundary conditions satisfied on the perturbed interface given by Eq.(l) are 



drj drj 

"777 ~t" VhxT{ 

at ox 



+ V hx— = Vhy (2) 



g^( V hx ~ V [x ) = V hy - V ty (3) 

where {vh,i) x , y are the velocity components of the heavier (density ph) and lighter (density pi) fluids 
occupying the regions y > and y < respectively. 



Assume single mode potential flow for the upper fluid governed by the potential 

h ( z> y } t) = a(t) cos{kx)e~ k{y ' r ' o(t)) (4) 

with « te = -^and v hy = -^ 

Since both V 2 Vh x = and V 2 ^ = for potential flow, the equation of motion of the upper 
incompressible fluid of constant density ph and coefficient of kinematic viscosity Vh (fih — PhPh) 



P h \~fo + (vh-V)vh] + V(ph + gp h y) + iihS 2 v h = (5) 
leads to the following integral for the irrotational motion. 

Ph[ ~ d ~w + 1^^ 2 + 9y] +Ph = m (6) 

For the lighter density fluid the flow inside the bubble is assumed rotational^ with vorticity Q = 
(-jfa - — ^q^-)z- The motion is described by the stream function ^f(x, y, t); 

vix = — — and vi v = — — (7) 
dy y dx K ' 

so that 



V 2 ^ = -u (8) 
Let x( x ,V,t) be a function such that 

V 2 x = -U3 (9) 
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Hence — x) is a harmonic function as V 2 (\I/ — x) = 0. Let $>(x,y,t) be its conjugate function 

a$ _ dx 

dx dy dy 

d± = _d^ + dx (1Q) 

dy dx dx 

i.e. that velocity components of the lighter fluid are 

Vlx 7{ 7{ T{ 

dy dx dy 
Vly dx dy dx 

The following choice is made for the stream function according to ^ 

*(x,y,t) = b {t)x+ [6i(t)e fc( ^" o) + u (t) / k 2 ] sin {kx) (12) 
with 

x(x,y,t) = u (t) sin (kx)/k 2 (13) 
Eq.(10) now gives 

$( x , y : t) = -b (t)y + h(t) cos (kx)e k{y - vo) (14) 

Substituting for the velocity components as determined by the velocity potential (f>h and the stream 
function \1/ in the kinematic boundary conditions Eq.(l) and Eq.(2) and equating coefficients of x 1 for 
i = and 2 we obtain the following four equations for nondimensionalized bubble height £i = krjo, 
nondimensionalized curvature £ 2 = r] 2 /k and b , b\ 



^f = -\(K* + 1 )& (16) 
kb 6^(26 -w) (17) 



V% (66 - 1) 

k 2 h (66 + 1)6-^ 
~ (66-1) 

Here 

k 2 a 



(18) 
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is the nondimensionalized velocity of the tip of the bubble, 



(19) 



t = t^kg (20) 
is the nondimensionalized time and 

(21) 

is the nondimensionalized vorticity. 

We need one more equation, i,e., the one for 6( r ) to complete the set of five necessary equations 
for the five functions describing the time evolution of bubble. This is obtained from the dynamical 
boundary condition as done below. 

Starting from the equation of motion of the lighter (lower) fluid with coefficient of kinematic 
viscosity v x (^ = p x v{) 

Pil-^j- + 2 V(^i 2 ) +oozx vt] + V{pi + pigy) - #V 2 U| = (22) 
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and using Eqs.(7)-(11) we get 

+ ^ gy] + - + + 

+ /"'O I -^' i!/) = P3) 
The last term on the L.H.S being the contribution from viscous drag. 

For rotational motion Bernoulli's equation does not exist except for special cases^. However for 
velocity potential and stream function as given by Eqs.(12)-(14) the integrations in Eq.(23) can be 
expressed in closed form^. 

Taking the difference of Eq.(6) and Eq.(23) at the interface of the two fluids and applying the dy- 
namical boundary condition, i.e, the pressure difference at the interface is balanced by the difference 
of viscous stress and surface tension^] we obtain the following: 

f ,du duo , . r , ~. . 
" J VA-Q-te - g^dy) = -[ph-pi] + f(t) 

where T is the surface tension and R is the radius of curvature at the interface. 

Substituting for fa, 1 i r , x, ®, using Layzer's approximation method^, expanding in power of x 
to 0(x l )(i < 2) and equating coefficient of x 2 we obtain after some laborious but straightforward 



manipulation the time dependence equation for 6( T )- The set of time evolution equations of the 
bubble are thus given by [together with Eq.(17) and Eq.(18)] the following three differential equations 
for £ 2 and 6 

f = - ^3(66 + D (26) 
^ = [-^2,r)-#- + 2(r - 1)(66 - 1)6(1 - l^fS)]" ' 



, r o> 2 - 5(66 + 1)^3 , ^ 1 
+ [ 71 7^ + w l 



(1-66) J ^(6,r) 
where 

iV(6, r) = 36(1 - r)6 2 + 12(4 + r)6 + (7 - r) (28) 



L>(6, r) = 12(1 - r)6 2 + 4(1 - r)6 + (r + 1) (29) 



,2 (ph-pi)g pi Hh ph v h k 

K = = ; s = — ; v h = — ; r = — ; c h = —= (30) 

T fi h p h pi y/kg 
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Ill NUMERICAL RESULTS AND DISCUSSIONS 

The time evolution of bubble is to be worked out by numerical integration of Eqs.(25)-(27) by 
employing Runge-Kutta-Fehlberg scheme. To this end it is necessary to know the dependence of the 
vorticity Cj(t) on r. We choose Cj(t) in the following form so as to have a time dependence having 
close resemblance with the simulation results^. 

"0") = 1 I x [tanh(r )(l + tanh(r)) + tanh(r - r )) (31) 

1 + 2 tanh(ro) 

Note that Cj{t) increases from £>(0) = and gradually builds up to a saturation value Cj c . The 
constants r and Cj c are adjusted so that the time dependence has approximate qualitative agrement 
with simulation results given by Fig.4 of Ref.^. This may be seen from the graph of Cj(t) plotted 
against r in Fig.l. 

In Fig. 2 are shown the nondimensionalized height (£i), curvature (£2) and velocity (£3) of the tip 
of the bubble and plotted against time r = \fkg. 



(a) . The first two terms [-Nfa, r)j^ + 2(r - 1)(6& - 1)6(1 - 12£f|)]-^ give the classical 
bubble velocity based on Goncharov model^ modified by effect of surface tension. 

(b) . The third term [ * "(fjg^ 3 + uj] ^ is the contribution from vorticity accumulation 
inside the bubble^ causing enhancement of the bubble velocity. 

(c) . The fourth term — 2cj l r[(s + 1)(1 — 12£f) + 4£ 2 (s — 1)] D $? ^ accounts for the lowering of the 
bubble velocity due to viscous drag exerted by the denser fluid. 

(d) . The last (fifth) term 2cf l rsu(l + 2£ 2 )(1 — 3£ 2 )/-D(£2 ; r ) gives the vorticity dependent contri- 
bution that diminishes the friction drag; such a possibility was suggested by Ramaprabhu et al.^D 
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Eq.(26) and Eq.(27) show that as r — > oo the bubble velocity reaches the following (nondimen- 
sionalized) asymptotic steady value with £ 2 — > — I (A— Atwood number), Cj{t — > oo) = Cj c and 



It is be noted that this result agrees with the saturation value of Betti and Sanz^J when viscosity 
is neglected(c/ l = 0) and with that of Sung-Ik-Sohn- - when vorticity is neglected. 

In Fig.2 the physical process(a) is represented by plot I. It is the Goncharov (classical) model 
result 



In absence of viscosity the vorticity aided time development of the bubble velocity is given by 
plot II in Fig.2. This represents the combined effect (a) + (b). As Cj(t) increases gradually from 
u)(t = 0) = 0, plot II indicates that so does the bubble velocity £3(7") from the initial value £ 3 (r = 0) 
remaining greater than the [^{r)} c i ass i C ai but close to it as long as the vorticity does not increase 
appreciably. This continues as r increases till Cj{t) builds up sufficiently as indicated in Fig.l with a 
consequent rapid growth of £3(7") above the classical saturation value [^{T)] c i assica i as given in plot I 
of Fig.2. As r 00 the vorticity aided bubble velocity approches the asymptotic value 



Cj{t — > 00) = 




(32) 




(33) 




(34) 
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As the bubble rises the viscous resistance of the upper (heavier) fluid depresses the vorticity aided 
bubble velocity. The magnitude of this depression depends on Ch = Ph^h ■ This is shown by plot 
III in Fig. 2 representing the cumulative effect (a) + (b) + (c). During the initial stage, Plot III is very 
close to plot I ([£,3(T)] c iassicai) due to vorticity and depressing effect of viscosity. But as Cj{t) becomes 
sufficiently large as r increases plot III shoots above plot I and attain the asymptotic value 



However due to accumulation of vorticity inside the bubble the viscous drag is reduced. This is 
indicated by the presence of the term sChCj c . The possibility of such an eventuality was suggested by 
Ramaprabhu et. al. and this is exhibited by the last plot IV in Fig. 2 representing the cumulative 
effect (a) + (b) + (c) + (d) which lies above plot III but below plot II. Thus the following inequalities 
hold among the asymptotic values: 

( V b )asymp ^ { v b)asymp < { V b) a symp (36) 

The leftmost inequality shows that the viscous drag exerted by the denser fluid as the bubble pene- 
trates through it is reduced by the viscosity effect on the motion of the lighter fluid due to vorticity 
accumulation inside the bubble. 
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Figure 1: Vorticity Cj{t) plotted against r with saturation value u c =2 and parameter r =8. 
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Figure 2: r = 5, s = 1/5, u c = 2, initial value £i = —£2 = £3 = 0.1. Black line 
(cD(r) = 0, Cft = 0), blue line (cD(r) in Fig.l, Ch = 0), dotted line (u>(r) in Fig.l, 
0^ = 0.2) and dash line (cD(r) in Fig.l, = 0.2) represents plot I, plot II, plot III and 
plot IV respectively. ^ 



